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STABILITY  OP  A COMPRESSED  ROD 
Karl  Kreutzer 

Introduction 

The  elasticity  problem  of  a break  £0  is  dealt  with  here 
from  the  viewpoint  which  considers  that  a break  process  is 
inherently  a problem  of  elasticity  theory  of  finite  deformations. 

It  is  assumed  that  the  tensions  which  appear  are  within 
the  validity  field  of  the  expanded  Hook  law  for  finite  deforma- 
tions. Included  here  is  the  limitation  of  fixed  mass  and  surface 
forces,  that  is  independent  from  deformation. 

The  method  used  is  an  energetic  one  which  takes  into 
account  the  finiteness  of  deformations  along  all  three  coordi- 
nates, as  developed  by  E.Trefftz  [2}  and  reported  at  the  inter- 
national Congress  for  Technical  Mechanics  in  Stockholm  (1930). 

The  following  nomenclature  is  used: 
yn  r..i  coordinates  of  points 

A,,  A',,  xt  components  of  mass  force  per  volume  unit  of  the 
non-deformed  body  along  the  three  directions 


components  of  surface  force  per  area  unit  of  the 
non-deformed  surface  along  the  three  directions 

m " displacement  components  in  the  three  directions 


components  of  state  changes  from  equilibrium 
(of  disturbances) 


I 


* 


E internal  energy  of  the  whole  body. 

An  elastic  state  is  a stable  equilibrium  state  when  with 
each  finite  displacement,  compatible  with  the  geometric  condi- 
tions the  increase  in  internal  energy  is  larger  than  the  work 
available  from  external  forces.  Thus  when 

1 A’>  Iff  v Xy  ds'"  dx“'  dx"‘  + (f  1 1 y Ah'"  do. 

* « 

Developing  both  sides  according  to  exponents  of  6„,r*  and  their 
derivations  gives: 

JE-  = tn+  &’E  \\  XIrAutrtdo  . . (1|. 

* » 

On  the  right  appear  no  exponents  of«»H‘”due  to  the  limitation  to 
fixed  external  forces  A,  and  I,  • 

Should  now  under  any  conditions  compatible  kinematically 
with  a,,"’*  the  left  term  be  larger  than  the  right  one,  then 
the  linear  terms  must  disappear.  Thus  is  should  be 

«£=■-  fff  -A",  6 n ,r*  dx'"  dx'*  dx*  -f  (f  X IyAo,r‘  do. 

9 # 

The  content  of  this  equation  is  described  as  the  "principle  of 
virtual  displacements".  It  expresses  that  for  each  virtual 
displacement  from  equilibrium  the  change  of  internal  energy 
equals  the  work  of  external  forces. 

Should  the  equilibrium  state  be  stable,  then  the  quadratic 
members  of  the  left  side  of  (l)  should  outbalance  the  members 
on  the  right  side  of  (1);  as  a consequence  of  limiting  ourselves 
to  constant  external  forces,  the  work  of  these  forces  is  spent 
by  the  linear  members.  The  stability  condition  is  thus  reduced 
to  **  0 *1 
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The  stability  limit  is  reached  when  for  at  least  one 
system  of  displacements  S u,r‘  the  second  variation  of  the  internal 
energy  disappears,  that  is  3*B-n  • 

These  very  general  points  will  be,  in  what  follows, 
applied  to  the  stability  problem  of  the  elasticity  theory; 
particularly  to  the  case  of  obtaining  break  loads  of  a compre- 
ssed rod. 

I.  Elasticity  theory  of  finite  deformations 
l.The  distorted  state 

To  define  the  particles  of  an  elastic  body  "substance 
coordinates"  are  used.  These  are  denoted  in  a non-deformed  state 
as  rectangular  normal  coordinates  and  indicated  by  In 

a distorted  state  they  become  curved  coordinates  correspon- 
ding to  each  mass  particle. 

Any  mass  particle  before  distortion  has  the  coordinates 
*"’•  *"’•  •**"•  It  includes  a point  P towards  which,  in  fixed  space 
intersection  lead  three  vertical  vectors  0,.  tt,.  C,  . and  the 
location  vector  . A neighboring  particle 

with  coordinatesar<,,+-«/.r"V"  + dx1”  *“*  + */*'**  includes  a point  Q. 

If  is  a vector  of  P according  to  Q then  rfi  = and 

the  linear  element  is 

• «• 

For  the  coefficients  of  the  linear  element 

under  normal  coordinates  the  following  scheme  applies: 


I 


ill 

0 

0 


ii 


0 

0 

l 


Through  an  elastic  displacement  & JV”  a,  the 

body  is  distorted.  The  point  P u)  takes  a new  position 

and  its  neighboring  point  Q a new  position 

(?(*  + rft  — i + «lj  M»  + «!*>)  . Besides, 

between  the  components  £"’»  of  the  location  vector  t 
after  deformation  and  the  components  *,r* of  the  location  vector  f 
before  deformation  there  is  a relationship 

= X,r'  + u’9*,  d(,r>  = d***»  + du «*\ 

The  linear  element  after  distortion  is 


and  where 


do,^£rr,l>d*‘"dxM, 

» p 

yd  ('*'  d _ d , d uM  , y di*'*’ 


> 


p 

The  linear  element  ds  has  thus  been  distorted  into  a 
linear  element  <*•*  • Comparing  both  changes 

G.„  for  the  linear  elements  coefficients,  we  get: 


ium  iwH  r 


The  axe  the" distortion  magnitudes"  describing  the 

tension  and  angular  changes  suffered  by  each  volume  element. 
Due  to  substitution  rules  there  are  only 

six  distinguishable  distortion  magnitudes.  These  form  a tensor; 

it  corresponds  to  each  point  of  on  elastic  body  and  is  its 
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symmetrical  distortion  tensor 


I 

' 


The  linear  element  do'  is  expressed  in  normal  rectangu- 
lar coordinates  as: 


1 

' 


where 


do'  ~ (1  + Y„)  dx“"  4 (1  + Y„)dx'n ’ 4 (1  + r„)  dx'"' 
f 2y„  dx dx'"  +l2y„  dx1"  dx'"  4 2y„  dx'"  dx'", 

- 2 0 x'" + Id  x'")  + Id  x'“)  U *"7 
d «"•  d «**'  , d n'"d  »*",  d «wd  «•"  , d «"'d  n,M 

/u  A .<«  i A .<»  ' .'l  r"1  A V’>  + d ^'"d  x***  rd  x^'d  x‘" 


< ■2*  • 


etc,  through  cyclic  substitution. 

These  non-linear  equations  transform  into  linear  equations 
of  the  classical  elasticity  theory  wnen  the  products  and  squares 

& If'** 

of  may  be  neglected  when  compared  to  linear  expressions. 

2.  Tension  state  and  the  equilibrium  condition 

In  order  to  describe  the  tension  state  of  a mass  particle 
f",  s*‘)  we  rega;r(}  the  square  angle  paralelopiped 

which  in  the  non-deformed  condition  is  formed  by  the  elements 
parallel  to  be  respective  axis  element^  d*'",ds'n,ds'" 
from  the  latter,  after  deformation  a general  paralelopiped  is 
formed  having  edges 

The  vectors  «/.=~^r„t,  which  give  the  direction  and  the 
irjrease  ratio  after  defoliation,  are  called  "grid  vectors". 

When  the  force  *"*  4*'" 4x‘"  acts  in  the  direction 

of  growing  <**"  boundary  surface  of  the  parallelopiped 
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i 
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, then  p'»  is  called  the  tension  vector  for  surfaces 

x^  = constant.  The  sane  corresponds  to  the  remaining  surfaces. 

,,r’  is  the  tension  vector  for  the  surface  element  •r,*',  = constant 
and  signifies  a force  per  unit  of  the  non-def ormed  surface. 

Each  of  the  three  tension  vectors  may  be  decomposed  according 


to  the  grid  vectors 


!m  - 11'  >•  Cm 
»/• 


;sAare 


by  which  nine  tension  componentsAare  created,  which  completely 
describe  the  tension  state. 

The  equilibrium  against  torsion  in  any  chosen  direction 

demands  the  disappearance  of  the  momentums'  sum  of  all  the  forces 

acting  upon  the  paralelopiped  considered  at  its  center,  thus: 

r.V  ^ l'(tr  X I'")  <lx"’  <tx”‘  ,tx,u  --o. 

t 

After  introducing  the  reciprocal  vectors 


of  grid  vectors  e’.e’.f*  according  to  formulas  C33 

and  keeping  in  mind  (3) 

£<»(*.»  fc*q  + c'”  (/."  - 1:'*)  -h  f'*‘  (J.”  - k”)  = o . 

. m •'  - « 

This  vectorial  equation  is  valid  only  when  ^ 

' !•  = * . 

Following  this  Cauchy  reciprocity  law  is  valid  for  l'* 

The  number  of  tension  magnitudes  is  reduced  from  nine  to  six. 

For  the  equilibrium  against  a displacement  in. .any  direction 
it  is  necessary  to  consider  the  action  of  force)  I'"  dx'"  ix'*  < 

(T ) * 

the  paralelopipeds'  surface  element  xv  /=constant  on  to  the 

t 

surface  element  •+■  ***"  — constant;  the  acting  force 


is  then  ?"**»*** 


♦J5  •"••• 


'lx,nflx"' 


The  force  excess  amounts  for  tiis  surface  pair  dj„,  dy  </*••'. 
Through  cyclical  substitution  follow  the  force  excesses  for  the 
remaining  surface  pairs.  A resulting  d'H  of  the  tension 
forces  is  obtained  which  is 

d Vr* 

</#  = - •!*'». 

Besides  a mass  force  acts  upon  the  paralelopiped 

dix  = ih/j-'"  <is-‘  <iy  - A / c, oy  dy  dy, 

i 

where  V is  the  mass  force  per  volume  unit  of  the  non- 
-deformed  body  and  its  component  in  the  direction  of  the 

th  grid  vector.  The  equilibrium  against  displacement 

requires 


<f:K  ( d M «>, 


from  here  after  abbreviating  through  the  product  «/.*-•»•  dy 


\ • «'  l*" 
T ’ 


o. 


Figure  1 
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Using  formula  (3)  and  differentiating  it  follows: 


■finis  it  is 


l/‘  > | 

where  I U I 

At  the  same  time: 


' ’ v ’ ^ . \ * \ ■ , t'  f <i  » , 

T - ! ' "a..."  I — o. 

• ft  1 •*  , 

•'  »<<  \ ' («  I ) 

■ tV1'1  T I A / r*  ’ 

are  Christoffer1  s symbols  of  second  class. 

\t  fit’’"  .«  , /«»i  n ) 

T T '"‘,Cr  V0 


Should  the  base  vector  be  everywhere  called  r*  , aniindex  ex- 
change has  to  be  effected.  If  we  decompose  this  vectorial 
equation  into  its  components  we  get 


T^'+T  r‘"{f)  +r'u-° 


(41. 


This  equation  delivers  for  each  of  the  three  directions 
h= 1 , 2 , 3 a partial  differential  equation.  The  middle  member 
considers  the  bending  of  the  coordinate  curves,  however  the 
diferentiation  in  (4)  follows  the  substance  coordinates  which 
are  curved  in  the  deformed  state. 

3.  Internal  energy 

The  approach  up  to  now  gives  each  point  a distortion  and 
tension  state.  To  represent  the  interconnection  of  both  we  have 
to  consider  the  internal  energy.  Its  existence  follows  for  all 
reversible  static  phenomena  from  thermodynamic  considerations. 
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The  following  meanings  are  used: 

*=*I«i*»f.lr=|/|/V*j  volume  of  each  grid  unit 
dl  i,ix  voiume  0f  infinitesimal  paralelopiped 

e f V o/am  drr/t  - /ion  a/s^ormc/t// 

n internal  energy  of  each  grid  unit  -600*  o 

tPK  - eilV  internal  energy  of  volume  dV 


In  order  to  represent  the  internal  energy  e as  a function 
of  the  distortion  magnitudes  y,,.  , we  have  to  find  the  in- 
variants of  the  distorted  state.  An  absolute  invariant  for  all 


values  of  paramater 


A i. 


\G.h jly’ " 

By  developing  the  determinants  by  exponents  of A a Jd 
degree  function  in  results.  Since  this  is  an  invariant 
for  all  y\,  the  coefficients  of  the  cubic  form  must  also  be 


for  all the  coefficients  of  the  cubic  form  must  also  be 
invariants.  Thus  three  invariants  are  formed: 

1 

' (-{n’)+y„ tdMo„  -a „9)d-/a4(6llr;..}  «„*) 

0,  ij. 

7* ■ (i,  r l • y***) -1- - ;>,.»)  + «„ (y, , y„ - y„») 

+ + - d3l(;'a,y„  y,,)} , 

/ !;'•  "! 

* j r;, 1 * 


i 


I 


For  I]_is  linear,  Ip  is  quadratic  and  I7  cubical, 

further,  Ip  follows  from  1^  through  a substitution  of  G with  ~Y 


For  rectangular  normal  coordinates  it  is: 

/ ••  i v ' ■■ 

‘ III  T-  y-y.  r , J,, 

h r„’/„  \ yn?» y..,‘ 


What  is  to  be  now  attached  for  e ? From  numerous  possibilities 
we  take  the  simplest  which  corresponds  to  the  classical  theory. 
Thus^* £ is  represented  by  the  simplest  quadratic  invariant 


of  the  distorted  state: 


1 '••"  I"  , , U I I 

1 i/v,  12  '• 


The  forward  factor  is  as  a volume  ratio  of  the  paralelopiped 
between  the  non-deformed  and  deformed  states  - an  invariant. 
The  constants  Of  and  A are  the  two  independent  elasticity 


constants  of  the  classical  theory 


(i  tn  1 a 

a=z  *,„  2’  P'~ 


G means  here  the  thrust  modulus  and  m PoXisson's  number  of  the 
cross-sectional  contraction  of  the  material. 

The  internal  energy  for  each  grid  unit  is  thus 

l*~  pljj (5). 


From  the  condition  that  e has  to  be  positive  it  follows  just  as 
in  the  classical  elasticity  theory 

» 

For  this  span  (5)  is  applicable,  since  («*)  • is  a positive, 
definite,  quadratic  homogeneous  form  in  r*r  • 
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For  rectangular  normal  coordinates  it  is 


2 t i'»»t  fllVn y»t  H'sjj'iitj’uj’ii  )*i*f  7»i*  y,,’)- 

The  infinitely  small  deformations  may  easily  be  arranged  into 
(5),  when  the  linear  expressions  are  introduced  for  the  distortion 
magnitudes. 

4.  Tension-expansion  equations  and  Hooke's  law  extended  to 
finite  deformations 

The  equilibrium  conditions  (4)  are  not  sufficient  alone 
to  determine  all  tension  and  distortion  magnitudes.  For  this  we 
also  need  the  relationships  between  the  acting  forces  and  the 
deformations  caused  by  them. 


>0.  


/ 

'A+Kar">Jd,«,- 

0.7. 


Figure  2 

As  in  the  classical  theory  a possibility  is  available  when 
observing  the  internal  energy,  to  deduce  tension-expansion  equa- 
tions, if  we  part  from  the  internal  energy  of  each  grid  unit. 

We  will  take  a real  displacement!^  and  add  to  it  an  addi- 
tional displacement  Sv  and  the  Yi>/<  will  change  by 
The  increase  in  internal  energy  is  equal  to  the  work  exerted 
by  the  tension  forces  upon  a volume  element. 
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Point  A suffers  an  additional  displacement  <»-/•'■•  -j 

and  point  B '**’  ' f"  . With  this  additional  displacement 

A ♦ L| 

an  internal  energy  change  is  obtained  f'"-  >”„.*/••■ '">h 
Through  cyclical  substitutions  internal  energy  changes  are 
obtained  for  the  remaining  two  directions.  Their  sum  comprises 
per  grid  unit: 


v s', 

r t'  r"' 4 

* if  • • 


now  «>n 

•»  L<-  * 


/•...  A/v( 


and  because 


it  is  Ar  ■ 1 


The  following  representations  are  obtained  for  the  internal  energy 
change : 

,V,i  r*'\v‘  : I 

On  the  other  hand  based  on  formula  (5) 

('tic.)  '(>(»•'').  ,<■(").  <'('•').  ... 

...  !*  ,v  *':*»* 4 *y  • t*«  • 

<’  Z»f  *'/1«  * * If  1 / #■»  1 

Comparing  both  expressions  we  get  the  six  equations 


b" 

•V,.  ’ 

* 1 /’it 

•>,'  (<!  r) 

»';(*  r | 

b” 

i'r,t  ' 

/.SJ 

" »r. 

b>" 

2 I'M 
'Vh  ” 

<-  ■ «'•.*» 
^ )’«l 

These  are  the  tension-expansion  equations  of  the  finite  defor- 
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mations  theory.  They  express  that  the  tension  magnitude; 
may  be  gotten  through  a partial  differentiation  of  specific 
internal  energy  fe  &)  according  to  the  distortion  magnitudes 

. They  are  valid  in  form  (7)  even  when  not  starting 
from  normal  coordinates. 

According  to  (5)  the  internal  energy  (eej  is  a homogeneous 
quadratic  form  of  the  distortion  magnitudes  . The  A 

are  according  to  (7)  linear  homogeneous  functions  of  the  distor- 
tion magnitudes  and  are: 


“ [u  [x  (o„ (i:u  a,,1)  i-('nyn)], 

V I l,rl‘  1 

T.  ,['i  \ | 

1 In  /i  1 


. . (8) 


etc,  through  cyclical  substitutions. 

These  equations  are  Hook's  law  extended  to  finite  deformations. 
For  rectangular  normal  coordinates  it  is 


k"  — <i\ 

-U 

<i>  \ 

k"  — Gy, 

1 / i i ' 

m 2/’ 

k”  — 0 I 

(;'«  + 

\ 

m -2/ 

. k” -«>•: 

k“  ( ; 

</>  | 
m 21’ 

A” -Gy. 

where  *»r.,+r,.+  r « 

Equations  (9)  convert  into  the  classical  tension— expansion 

equations  when  infinitely  small  deformations  are  allowed. 
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Then  with 


a 

«.r, 

D H 
i)  x' 

(3) 

[ ;n  w? 

'-'ft 

k 

13 

<•  ("//-! 

'’.rb 

);■ 

r;t  . 

<'<  (<-z  ! - 

/ (->  \ 

k 

ii  . . 

<i(Wy  1 

'O- 

Equations  (9)  supply  the  connection  between  the  tension  components 
and  the  distortion  magnitudes.  If  the  tension  component's 

values  from  (9)  are  put  into  the  equilibrium  conditions  (4) 
we  obtain  for  the  determination  of  the  equilibrium  the  necessary 
differential  equations  for  the  displacement  components  ; j 
. For  these  differential  equations  (exactly  as  in  the 
classical  elasticity  theory)  correspond  limiting  conditions 
which  ma  refer  to  displacements  or  to  surface  forces. 

5.  Stability  of  equilibrium 

Above  we  presented  equations  which  serve  to  ascertain  he 
equilibrium.  For  the  study  of  stability  it  is  anticipated  that 
these  equations  be  integrated  so  that  we  may  know  the  tensions 
and  displacements  of  the  equilibrium  studied  for  stability. 

To  apply  the  theory  we  sha  1 satisfy  ourselves  with  a relatively 
simple  particular  case  - a compressed  rod-  whose  tension  and 
displacement  states  may  simply  be  overlooked. 
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To  deal  with  the  stability  we  have  to  compare  the  internal 
energy  of  tie  body  in  equilibrium  state  with  tie  internal  energy 


» 


-4 


in  a neighboring  state.  Thus,  besides  considering  the  equilibrium 

displacements  ^ ^ we  have  to  consider  the  neighboring 

| .«“»  f A h"\  m"*  + Am"',  m"’  + A m"’J 

displacements  amd  the  expression  for  the 

internal  energy  becomes 

E = ///  (c  f)  d xu'  tlx,n  dSn 

in  the  periphery  of  the  equilibrium  state  by  exponents  of  and 

eventually  developed  by  their  derivation. 

According  to  the  statements  in  the  introduction  the  equi- 
librium will  be  stable  when  for  each  allowable  ( thus  compatible 

x (yj 

with  the  geometric  conditions)  system  of  "disturbances" 
the  so  called  second  variation  (which  when  developed 

contains  the  quadratic  members)  is  always  positive. 

If  we  describe  the  differentiations  by  the  following 
indices,  for  exampel  «,*r*  = , then  the  second  variation  will 

acquire  the  following  form  4) 

* R~  w { ii  2 ' i<  ® (* + 2 fi  (« r;  ,r + 2^4%  $ + **  *i  A a 

+ V7 £ !2  /*(»  -t-l*’)'  + 2 ft ( ,4*1  + ,, ,)’  + 4 o (n<"y  + *»]  A «[“  A 4 

+ 2 [2(o  /?)(!  f-  4>)  + *•**'!  A m!/'1  <)4> 

'•  0 (* + + - /*(l + «!*')  4*  • + 2 /M*<1 + a,  , A)j  A a< 

+ 2',j4(«-  *>«  ,,«*•»!  A „!*'  A 

+ 2l4(a  /*!( I bH/iuf’ 

4 2’ 12(2 a -//)  «£>  ,,,4"' 

/»,/.,  VI 
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The  indices  go  through  values  from  one  to  three  and  two  differently 
numbered  indices  may  never  acquire  the  same  value.  Values  above 
three  of  the  indices  are  to  be  reduced  by  modulus  three. 


Nov/  the  matter  is  to  decide  in  special  cases  if  such  allow- 

v (a 

able  disturbances  0/1/  nay  be  found  for  which  the  second 
variations  becomes  negative  or  if  it  remains  positive.  Nov/  then 
the  displacements  and  tensions  of  the  eouilibrium  state  are  to  be 
regarded  as  given  and  we  should  now  look  for  the  "most  c angerous" 

f 

disturbances,  such  0 which  v/ould  possibly  allow  the  second 

variation  to  be  negative. 

6.  Determination  of  the  stability  limit 

To  determine  the  stability  limit  of  ah  elastic  equilibrium 
state  v/e  need  to  consider  the  second  variatinn  of  the  internal 
energy.  It  may  be  obtained  from  expression  (10)  if  we  introduce 
into  it  the  values  corresponding  to  this  special  case  for  dis- 
placement and  tension  components.  Since  the  quadratic  form  in 
SJ"  under  the  integral  and  its  derivations  may  always  be  written 
as  a difference  of  two  positively  defined  forms,  we  should  write: 

A'  #!k  I)  — if 

• ' I »• 

Where  Q-.  and  Q0  are  integrals  over  positively  definite  quadratic 

j (l>) 

forms  in  the  derivations  of  O W' 

To  find  if  the  difference  could  become  negative,  v/e 


write: 


~ V,  (i  — i) , 

v» 


The  smaller^  becomes  the  higher  the  danger  f or^jr  E to  become 
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negative.  The"most  dangerous"  disolacement  from  equilibrium  state 
is  such  for  which  /'i  becomes  a minimum. 


If  ^ becomes  a minimum  then 
must  disappear,  thus  it  should  be 

AQ.-iiQ, 


,•> •= 


b,Ao, 
o = 


III' 


for  all  allowable  variations  cJ^cf ) 

Equation  (11)  has  a solution which  may  be  larger, 
equal  or  smaller  than  one.  The  case^  — / gives  the  stability 
limit. 

Since  we  are  asking  only  for  a stabil  ity  limit,  the  pro- 
blem may  be  simplified  further.  Equation  (11)  says:  for  the 

r (o  {■  (v  f (i/ 

"most  dangerous"  displacement  O'L* ^ t U fr0m  equilibrium 

state  is  for  all  changes  $ ) 

Now,  instead  of  finding  the  value  A for  all  the  given  forces 

and  then  ask  for  which  forces  /[-  / , we  may  directly  put  /?-/ 

.Then  we  will  seek  the  forces  and  the  corresponding 

displacements  fid*,  j ^ d from  equilibrium  state;  for 

any  6C^^VJ)  it  will  be 

6Q,~*Q, 

This  equation  is  actually  Jacob's  criterion  for  the  access  of 

stability  change.  It  allows  the  obtainment  of  the  stability 
limit  for  any  equilibrium  state. 

It  should  be  noted  that  the  stability  limijr  may  also  be 
determined  from  the  isoperimetric  problem:  among  all  the  allow- 
able displacement  variations  from  the  equilibrium 

state  the  "most  dangerous"  are  those  for  which  the  integral  Q* 
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becomes  a minimum  for  a neighboring  condition  Q~=l.  With  **- 
as  the  Lagrange  factor  for  the  neighboring  condition  it 
becomes  identical  to  equation 

6Q  =A<5Q 

which  agrees  with  (11). 

The  stability  limit  criterion  expressed  here  in  connection 
with  the  variation  calculation  may  also  be  converted  into  a 
homogeneous  system  of  differential  equations  and  homogeneous 
limiting  conditions  for  the  "most  dangerous"  variation.  Such  a 
homogeneous  problem  has  several  solutions  different  from  zero 
when  the  parameters  characterizing  the  equilibrium  state  (loads) 
acquire  definite  critical  values  (break  values). 


rod. 


This  method  will  be  used  below  with  the  example  of  a compressed 


II.  Stability  of  a Compressed  Rod 


7.  Equilibrium  state 

To  carry  out  the  Trefftz  method  in  order  to  find  the  stabi- 
lity limit  of  a compressed  rod  the  definitions  should  be  fitted 
to  this  particular  problem. 

The  coordinates  are  designated  x,y,z.  The  coordinate  system 
is  layed  out  so  that  the  z-axis  coincides  with  the  rod  axis  and 
the  remaining  axis  fall  within  the  cross-section  of  the  rod. 

The  equilibrium  displacements  are  indicated  with  capital 
letters  U,V,3ftC  Variations  of  U,V,\8(will  be  indicated  by  small 

\jj 

u,v,*  and  the  variations  of  these  latter  will  be  indicated  by 

Jw  S'O'  . S L*r 

/ / 

The  rod  considered  should  be  "fixed".  This  means  that  all 
points  of  the  lower  and  all  points  of  the  uppper  cross-sections 
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have  suffered  the  same  vertical  displacement;  thus  V should  be 
predetermined  for  the  cross-sections  of  the  extremes.  From  all 
the  disturbances  of  the  equilib  ium  state  only  those  are  acceptable 
whose  w=0  at  the  frontal  surfaces  of  the  rod,  and  thus  also 

The  rod  is  held  by  its  lower  end  (z=0)  and  compressed  by  a 
portion  AL  by  a force  P=p*f  (f  is  the  cross-section). 


then  the  displacements  are 


U’—  u:r, 


and  these  may  be  neglected. 


The  connection  between  the  expansions  ax,a^.,az  and  pressure 
p is  obtained  from  the  tension-expansion  equations  (9)  which  now 


apper  as 


a--*-'-  a i 


<t>  \ 
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^ -3J  = o. 
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Figure  3 
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The  equilibrium  conditions  (4)  are  (as  in  any  tension  state) 
self-fulfilling. 

8.  Second  variation  of  the  internal  energy 

The  homogeneous  equilibrium  loses  above  a "critical’magni- 
tude  of  pressure  p its  stability. 

To  find  the  stability  limit  the  second  internal  energy  va- 
riation has  to  be  considered.  It  results  in  : 


/',  - |f|  {:2ii|(1  a.,r  «v!  4-(14-n.»/)'  «V4-.0  4 "4J  '.-'I 
-i-4ld  /»)  |(  I 4-0//I  «>  i’.i/  III  f-  it,,)  |l  H- if,  4(1  ^U;)|l4-«t)>i'.  n.r 

I /-*  |(l  -4  nxy  («.**  4-  «»’)+( 1 5 r’ 4-  «'-•')  + (l  -!■ «.!'  («’*’  4 iv 

4-  >|ii  . j.  !•„)  v„  ; :?  *,  I ! i -«l  I n.l(l  I 

\ l \ \ k«  ' 4"  z »'\  v)  d <4>'  ■ / ? . 


introducing  instead  of  disturbances  u,v,w 


i‘  1 1 i~ «,  j a . 
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and  with  a i • 

j-  ,T'S 


-20- 


, . I 11  l1'  III  I . _ _ 

•*  '•  \\\  | .1  M .((".«  i »■-/  i-  - a 1 1 «,  <’»  t i »•„  i< , 

I tv.  II y (i<„  + rxy  (»•.  (iT'r •}•  i<;)?|J  il  vil nils 


H;  I'/  , M*,’ 

(I  I-  ii v) 7 (I  ii «r  ^ ( I -4- 


>/.r  ill)  il z . 

The  first  integral  is  equal  exactly  to  the  form  change  work  of  the 
classical  elasticity  theory.  It  shows  the  form  modification  work 


involved  in  small  displacements 


when  these  alone 


are  present.  The  second  integral  contains  the  acting  pressure  p. 

And  so  the  second  variation  of  the  internal  energy  of  a compressed 
rod  is  written  as  a difference  of  two  integrals  over  positively 
defined  quadratic  forms  in  derivations  of  ^ , 

thus:  <2.  -&-L. 

9.Jacobi's  equations 

After  developments  in  6,  we  reach  the  stability  limit  for  an 

elastic  equilibrium  state;  if  there  is  a "most  dangerous"  variation 

u,v,w  for  which  with  every  acceptable  variation  j cf 

UiS 

and  for  the  case  of  a compressed  rod 
and  Q~  have  the  above  values. 

The  satisfaction  of  equation  (12)  for  any  variations  of 

u,v,w  leads  to  a system  of  three  partial  differential  equations. 

It  represents  Jacobi's  equations  adapted  to  an  equilibrium  problem 

in  a form  of  a variation's  problem  (6). 

To  simplify  the  derivation  the  crosslines  above  ^ ^ ^ are 

again  omitted.  The  partial  derivatives  of  the  integrants  of 

along  the  break  magnitudes  H etc  are  ;’‘n 

formal  accordance  with  the  classical  theory  designated  by  the  ten- 


sions 


etc. 
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Since  Q1  has  the  shape  of  Hooks  form  change  work,  the  first 
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variation  of  Q1  results  as  when  "•  W = «,  + *•*  + «» 
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ion  may  be  possible  for  all  acceptable  variations 


only  when  the  integrants  in  the  space  integrals  are  mutually  equal. 
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It  follows: 
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These  are  Jacobi's  equations. 

Limiting  conditions  belong  to  the  above  and  follow  from  the 
surface  integrals. 

a)  Integrals  touched  only  by(7^P  over  the  shell  plane  have 
to  dissappear.  It  follows  from  this  that  on  the  shell  surface 
following  equations  apply: 


Ox  CHS  (II,  .••)  - ,,  c.os  ( II , If)  (I 

Tx  U COS  (»,  J ■)  Oy  COS  (ll.  If)  0 

TT  _.  COS  (V.  ./•)  h Ty  2 COS  (,/,  If)  : (I 


(cos  II,  S =--{)) 


b)  The  integrals  oS:^6p.  and  U U taken  over  the  front 

I ^ 2, 


They  indicate  that  a force-free  shell  plane  corresponds  to  the 
disnlacements  u,v,w. 

' to, 

surfaces  must  be  equal.  Since  according  to  7,  for  all  allowable 
disturbances  and  their  variations  at  both  front  surfaces  are 
w=0  and  -^Ofor  all  points,  also  wx=0  and  w=0,  it  follows  that  for 
the  front  surface = Q and  So  f or  ^ upper 
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front  surface  expression  only 


remains  and  for  any  chosen  Su, 

J2-0. 

oz 

The  same  applies  to  the  lower  front  surface 

llv=  0. 

bz 

With  the  condition  w = 0 we  have  for  the  limiting  conditions 
at  the  front  surfaces 
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3,=°  od,J''  m—0,  »«,  — 0. 
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Finally  we  haveYthe"most  dangerous  displacement"  u,v,w  is 
sufficient  for  the  homogeneous  differential  equations 
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with  homogeneous  limiting  or  boundary  conditions; 
a)  for  the  shell 
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b)  for  the  front  surfaces 

" " • <* , i\  a , 

and  the  tensions  according  to  the  classical  theory 

...A'"  H \ a 'a  ti'\ 

. s). 

Equations  (13)  give  a simple  meaning  which  gives  the  connec- 
tion of  the  extended  theory  with  the  elementary  approximation 
theory.  When  in  the  equilibrium  the  equations  of  the  classical 
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are  put  as  the  volume  forces  we  obtain  equations  (13)-  This  re- 
veals that  the  integration  of  these  equations  is  identical  to  the 
task  dealing  with  finding  the  equilibrium  profile  of  the  rod 
when,  at  a proper  p-value,  the  loads  per  each  volume  unit  are 
proportional  to  the  magnitudes  '^£2  y 

When  the  volume  load  over  the  cross-section  of  the  rod  is 

■the' 

integrated  we  obtain  loads  per  length  unit  of  -fete-  rod 
P (PM  iY,ir-  1‘  ( '/  .1  y — p 

J A"  t — ■ (,  +fix),  (U,  . \ )<l  h - {l  4. tlz!  ’ \/‘dt  ~ +aj;  (if  ’ 

where  are  the  average  values.  If  these  average  values 

fl»  ? ' 

are  replaced  by  the  differential  ratios  ^ 

taken  for  the  rod's  centreline,  and  the  small  magnitude  ct  in  the 
denominator  is  neglected,  the  task  arises  to  find  the  equilibrium 
profile  of  the  rod  which  is  loaded  vertically  with  forces 


<Pi> 

> 

fl  9J 


t <!*  /> 
fl  zJ 


. This  is  exactly  the  content  of  equation 


hlwiv—  pwii  which  follows  from  the  ele- 
mentary theory. 

III.  Numerical  Results 

The  Jacobi's  equations  together  with  the  boundary  conditions 
of  page  361  (23)  tell  us  that  there  exists  a neighboring  equi- 
librium state  next  to  the  initial  conditions.  The  displacements 
that  transfer  the  initial  conditions  into 

the  neighboring  conditions  are  the  similar  factors//*  ^ etc. 

solutions  u,v,H'.  Since  the  Jacobi's  equations  and  the  boundary 
conditions  are  homogeneous  we  have  a proper  value  problem.  The 
proper  value  is  the" critical"  load  p. 

• -26-  V- 


In  this  section  we  solve  the  proper  value  problem  for  cir- 
cular and  rectan pillar  cross-sections.  In  the  fcrst  case  - by 
integrating  by  development  in  series;  in  the  second  by  the  Ritz 
method. 


10  Circular  cross-section 


Equation  (13) 
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is  transformed  for  cylindrical  coordinates.  This  moans  we  seek 
the  equilibrium  profile  of  a rod  loaded  per  volume  unit  by  forces: 


''  VTF  A- 


where  the  very  small  magnitudes  jOf  G g are  negligible  when 
compared  to  1 . 

When  the  cylindrical  coordinates  are  introduced  f "0  2 and 
the  displacements  are:  radial  £ , tangential  t and  ?27 

axial  4^,  the  problem  is  to  find  the  equilibrium  profile  for  a rod 
loaded  in  the  three  directions  by  volume  forces  ' p » t’ 

when 


are  the  cylindrical  coordinates  of  the  tension  forces. 
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The  equilibrium  conditions are 
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and  the  tension-eapansion  equations  are 
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By  eliminating  tensions  in  the  equilibrium  conditions  through 
a tension  expansion  equation  we  have  Jacobi's  equations  with 
cylindrical  coordinates 
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To  integrate  (14)  we  chose 
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It  is  sufficient  for  boundary  conditions  at  the  front  planes 
when  we  take 


For  functions^f/^^^/^/^^ythree  common  differential  equations 
follow  from  (14) 
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Where  the  vertical  lines  signify  derivatives  according  to 
and  q equals  p/G.  ^ 


Through  the  inotroduction  of  cylindrical  coordinates  it  was 
succeeded  to  reduce  the  partial  system  (14)  to  a system  with  three 
common  differential  equations  of  second  order. 

Equations  (17)  supply  the  cane  of  a break,  that  is  a linear 
bending  of  a solid  cylinder  and  a hollow  wall  column  for  n=l. 

The  equation  s are,  when  the  cross-contraction 

value  is  10/3: 
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The  boundary  conditions  for  the  shell  are  simplified  by  using 
(13)  in  equations 


7 r„P'  +9(P+Q)+:\yr„  /?-■*- 0,  | 

P+  <1  r„</—0,  I . • fill). 

1C  Vl'r=t)  ) 

and  in  case  of  a hollow  cylinder  there  are  three  analogous  equa- 
tions with  instead  of  7c)  . 

A complete  integration  of  (18)  leads  to  six  integration 
'constants.  The  correspond  with  six  boundary  conditions  (three  at 
the  inside  and  three  at  the  outside  boundaries).  For  the  solid 
cylinder  there  are  three  conditions  for  the  outside  boundary. 

For  the  three  missing  conditions  we  hove  that  for  tensions 
/ for  r-0  are  finite.  Functions  P,Q,R  fcr  YzO  c»***tt 


behave  in  a regular  manner.  This  regularity  requirement  supplies 
three  more  data  for  the  determination  of  constants. 

Equations  (18)  are  integrated  normally  by  developing 
exponential  series  for  the  solid  cylinder.  The  succesive  calcu- 
lations of  coefficients  delivers,  designating  the 

coefficients  c^,c0,c^,  the  series: 
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Further  calculation  is  possible  only  numerically.  Because  of  the 
complex  structure  of  (21)  with  respect  to  a,  load  values  are  used 
to  solve  the  series.  Since  q=p/G  and  the  push  module  G is  very 
high,  q must  be  small.  With  G=800,000  kg/cn^  we  obtain  for  the 

p 

loads  p,  which  are  within  Hook's  law,  p=800kg/cm  and  p=l,600kg/cm 
the  values  q=l/l,000  or  q=l/500. 

The  results  of  equation  (21)  were  compared  with  the  elemen- 
tary theory  where  the  break  lenghtk'^Ts  given  by 


The  result  is  that  there  occurs  only  a small  deviation  from 
Euler's  formula.  It  is 
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Thus  within  the  calculation  accuracy  Euler's  break  formula  is 
confirmed. 

11.  Rectangular  cross-section 

The  rod  has  a lenghth  L.  The  rectangle's  dimensions  are 
2a  and  2b (a  >b).  The  position  and  the  tension  relationships  are 
the  same  as  before. 


The  stability  limit  is  dominated  by  equation 
and  Qp  h0ve  the  meaning  of  page  360  (23). 


where 


-33- 
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The  Ritz  method  is  used  to  solve  the  equation. 

For  displacements  u,v,w  we  chose 

they  reduce  the  spatial  problem  to  a plane  problem.  It  satisfies  the 
boundary  conditions  for  the  front  planes,  where  for  z=0  and 
z=L  it  has  to  be  w=0,  u?=0,  =0.  The  integration  of  and  Qg 

over  the  rod  lenghth  L gives 
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For  A,B,C  we  chose 
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Introduced  in  (23)  we  have 
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Integrals  and  Qp  became  the  quadratic  functions  of  the  18 
Efri gfnYHjalt  coefficients.  The  equations  for  the  stability  limit 
AV.  = A Q„ yields  18  linear  equations  for  the  coefficients  according 


to 


. Thq  18  equations  system  is  distributed  into 


4 groups  v/hich  contain  the  coefficients  l- <*«;  A:  ftt 
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The  individual  groups  are  systems  of  homogeneous  linear  equa- 
tions where  coefficients  have  a non-zero  solution  only  when  the 
determinant  disappears.  If  we  take  out  one  group  and  equate  its 
deerminant  to  zero,  this  determines  the  break  load,  we  will  ob- 
tain the  coefficients  of  this  group  up  to  a common  factor. 

The  coefficients  of  the  other  group  become  zero  because  the 
determinants  of  the  other  group  for  the  break  load  value  found  do 

not  disappear.  We  are  interested  here  in  the  determinant  equation 
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which  supplies  the  break  load,  thus  the  question  arises:  which  of 
the  four  groups  should  we  take?  Naturally  such  for  which  the  lowest 


break  load  results; 


which  one  it  is  - is  not  difficult  to 


recognize.  If  the  rod  breaks  in  the  x direction  then  all  the 
cross-section  points  have  about  the  same  displacement  in  the 
x direction;  u cannot  be  thus  an  uneven  function  of  x.  Nov/  the 
only  group  which  has  no  uneven  function  for  u is  the  third  group. 
Putting  the  determinant  of  this  group  as  equal  to  zero  we  should 
obtain  the  break  load  value  which  causes  a break  in  the  x-direction 
The  fourth  group,  which  follows  from  the  third  group  by  substitu- 
ting x and  y,  yields  the  break  load  for  a y break.  Thendeterminant 
equation  is  the  same  as  for  the  third  group  when  a and  b are 
substituted.  Groups  1 and  2,  which  v/ith  the  substitution  become 
one,  give  a kind  of  squeeze  boundary,  but  this  value  has  no 
physical  meaning  because  these  phenomena  take  place  outside  the 
limits  of  the  expanded  Hook's  lav/. 

We  write  the  third  group  of  equations  as  follows 
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Since  the  coefficients  of  the  exponents  of  q change  with  its  sign, 
the  equation  may  only  have  positive  roots. 

For  a numerical  calculation  we  note  that  q as  well  as  k are 
very  small.  If  we  consider  the  limit  case  of  very  thin  rods 

O 

(we  let  k tend  to  zero)  then  after  multiplying  by  k we  are  left 
with  the  last  members  in  the  linear  and  absolute  members  which 
yield  the  equation  for  q 

13k2-15q=0 


From  this  - the  break  tension 
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Introducing 


r i’  10  . 

+ m T»- 


I «’ 

F ' 3 ’ 


i’k  =;»*  • f, 


we  obtain  exactly  Euler's  formula  for  the  tensioned  rod: 
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This  formula  is  thus  confirmed  for  tie  limit  case  of  very  thin 
rods.  The  lateral  relationship  n does  not  appear  explicitly. 

Strictly  speaking  there  is  a need  for  a n-depending  correction 

but  in  practically  important  cases  this  correction  does  not 
reach  any  noticeable  effect. 

I 

I . Summary 

As  an  example  of  the  stability  theory  for  elastic  equili- 
brium developed  by  E.Trefftz  the  present  work  calculates  the 
tv.  break  limit  for  circular  and  rectangular  rods.  The  first 
part  outlines  the  elasticity  theory  for  finite  deformations. 

In  particular,  it  is  shown  that  the  expression  f r the  elastic 
potential  (internal  energy  of  a volume  unit)  may  be  taker!  from 
the  elasticity  theory  of  small  deformations  when,  instead  of 
linearized  distortion  magnitudes  of  the  classical  theory  the 
real  distortion  magnitudes  for  finite  deformations  are  taken. 

(because  of  this  the  coefficients  of  the  linear  elements  change) 

From  the  obtained  integraal  for  the  overall  inner  elastic 
energy  this  work  develops  an  expression  for  the  second  variation 
of  internal  energy,  whose  sign  decides  the  stability.  Then 
following  common  methods  for  variations'  calculation,  the 
stability  limit  is  determined,  that  is  the  load  limit  above 
which  the  second  variation  is  capable  to  accept  negative  values. 

The  third  section  contains  results  of  numerical  calculations 
. . from  which  it  follows  that  within  the  limits  of  calculation 

accuracy  (up  to  fractions  of  one  percent)  Euler's  formula  is 
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bonfirned,  which  formula  ir-  obtained  according  to  elementary 
beam-bendinr  theory. 

I wish  to  thank  here  Proffecsor  Dr.  Trefftz  for  his 
encouragement  in  my  work. 
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4)  This  expression  was  first  used  by  Trefftz  in  1930  and  was 
given  to  me  privately. 
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